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ABSTRACT 

Aims. We intend to spatially and spectrally resolve the [O i] emission region in two nearby Herbig stars. 

Methods. We present high-resolution (A/AA = 80,000) VLT/UVES echelle spectra of the [O t] 6300A line in the Herbig Ae/Be stars HD 97048 
and HD 100546. Apart from the spectral signature, also the spatial extent of the [Oi] emission region is investigated. For both stars, we have 
obtained spectra with the slit positioned at different position angles on the sky. 

Results. The [Oi] emission region of HD 100546 appears to be coinciding with the dust disk, its major axis located at 150±11° east of north. 
The SE part of the disk moves towards the observer, while the NW side is redshifted. The [O i] emission region rotates counterclockwise around 
the central star. For HD 97048, the position angle of the emission region is 160± 19° east of north, which is the first determination of this angle 
in the literature. The southern parts of the disk are blueshifted, the northern side moves away from us. Our data support the idea that a gap is 
present at 10 AU in the disk of HD 100546. Such a gap is likely planet-induced. We estimate the mass and orbital radius of this hypothetical 
companion responsible for this gap to be 20 Mj llpiter and 6.5 AU respectively. 

Conclusions. Based on temporal changes in the [O i] line profile, we conclude that inhomogeneities are present in the [O i] emission region of 
HD 100546. These “clumps” could be in resonance with the suggested companion, orbiting the central star in about 11 yr. If confirmed, these 
observations could point to the existence of an object straddling the line between giant planet and brown dwarf in a system as young as 10 
million years. 

Key words, circumstellar matter — stars: pre-main-sequence — planetary systems: protoplanetary disks — planetary systems: formation — 
stars: individual: HD 97048 — stars: individual: HD 100546 


1. Introduction 

HD 97048 and HD 100546 are two of the most famous Herbig 
Ae/Be stars. The shape of the mid-IR spectral energy distribu¬ 
tion (SEP) sug gests that their circumstellar disks have a flared 
geometry (tDominik et al.lE()03l) . For HD 100546, broad-band 
images —sensitive to scattered light— of the disk structure 
have been obtained in the optical and near-IR JPantin et alJ 
2000i lAugereau et alJl200ll iGradv et all 1200 ll) . Furthermore, 


the disk around HD 100546 has been spatially resolved at UV- 
to-mm wavelengths in the continu um as well as in spectral fea- 


tures emanating from the disk (e.g. Wilner et alji2003; Liu et al. 

2003 

van Boekel et al. 

12004; 

Leinert et al. 2004; Gradv et al. 

2005 

). The disk’s appearance and position on the sky are hence 


relatively well-established for this source: the inclination of the 
system is f^50° and the position angle of the major axis is ap¬ 
proximately 150° east of north (e.g. lLiu et al.ll2003l) . The disk 
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of HD 97048 on the other hand has not been observed in scat¬ 
te red light, although the m id-IR (around 10 pm) observations 
of lvan Boekel et al . ( 1200 4) show that the source is resolved in 
both the continuum and the PAH bands. 

Herbig Ae/Be stars are by definition emission line stars. 
In a significant fraction of these objects, also forbidden transi¬ 
tion lines are observed. In a previous paper JAcke et alJll2005l. 
hereafter AVD05), we have studied the [O i] lines at 6300 and 
6363A. We have shown that for HD 97048 and HD 100546, 
the spectral line profiles are in agreement with the theoretical 
profiles emanating from excited neutral oxygen atoms in the 
surface of a flared, rotating passive disk. We have demonstrated 
that these atoms cannot be thermally excited, and suggested the 
excitation may be due to the photodissociation of OH and H 2 O 
molecules in the surface layers of the disk by the UV radiation 
field of the central star. 

Since both sources are relatively nearby and the emis¬ 
sion region of the [Oi] 6300A line is expected to be spa¬ 
tially extended, we have tried to obtain spatially resolved high- 
resolution spectra of the targets around this line. We present the 
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data set and reduction method in the following Section. The 
analysis of the data and the confrontation with the model are 
described in Sects. 0and0 In Sect.[5]we discuss new evidence 
for the presence of a giant planet around HD 100546, based 
on the results of the analysis. The mass and orbital radius es¬ 
timates for this object (M v i ane t = 20 Mj up iter and R = 6.5 AU 
respectively) agree strikingly well with previ ous st udies con¬ 
cerning the gap i n HD 100546’s disk jRoiiwman -glai J lizQQ-S i 
(Grad y et al.l2005l) and the spiral-arm structure in the outer disk 
flOu^en"et alihooll) . The final conclusions are summarized in 
Sect. |6l 

2. Data set and reduction method 

We have obtained a number of high-resolution (T/AT =80,000) 
VLT/UVES 1 spectra for HD 97048 and HD 100546, with the 
slit positioned at different angles on the sky. For HD 100546 we 
have taken spectra with the slit aligned with (PA = 150° east of 
north) and perpendicular to (PA = 60° E of N) the known major 
axis of the system, as well as at two intermediate PAs (15° and 
105° E of N). For HD 97048, no information on the disk’s incli¬ 
nation or major axis is present in the literature. Therefore, we 
have chosen random, equally spaced PAs of 0°, 45°, 90° and 
135° E of N. The slit length is 12", while its width is 0.3". In 
Table [U the observation log is presented. 

UVES is an echelle spectrograph, which implies that the 
raw 2D CCD frames contain one spectral direction (along the 
different orders, perpendicular to the slit) and one spatial di¬ 
rection (perpendicular to the orders, along the slit). The small 
width of the slit ensures the high spectral resolution needed 
to spectrally resolve the [Oi] line. In the spatial direction the 
spectra have a resolution of 0.182"/pixel, or equivalently ~64 
pixels over the 12" slit length. In a standard reduction scheme, 
the spectrum of the target is retrieved by integrating the flux 
in the spatial direction. However, we intend to use both the 
spectral and spatial information contained in the echelle or¬ 
ders. Therefore, we have taken the pipeline-reduced 2D echelle 
frames as our starting point. Spectra obtained with the same slit 
PA were averaged such that the peak of the spatial profile in the 
continuum of both frames coincides. In the following, we only 
consider the small part of the UVES CCD around the 6300A 
line in the averaged 2D spectrum at each slit PA. 

[Oi] emission of the Earth’s atmosphere, the so-called air- 
glow, can easily be removed from the spectra. The scientific 
target only fills a fraction of the slit, while the airglow feature 
is equally strong in all spatial directions on the sky around the 
source and hence fills the entire slit homogeneously. We have 
determined the spectral profile of the airglow feature by aver¬ 
aging 22 spatial rows at the top and bottom edge of the echelle 
order (i.e. outside the region of the slit which is dominated by 
the target). This average is then subtracted from all (spatial) 
rows, which results in a 2D airglow-corrected frame. 

Instead of integrating over the spatial direction and achiev¬ 
ing the “standard” spectrum, we have determined a normalized 
spectrum at each sampled spatial position. First we have deter¬ 
mined the continuum on the blue and red side of the line. We 

1 http://www.eso.org/instruments/uves/ 


Table 1 . The data set. For both stars, spectra were obtained at 
four slit PAs. 


Object 

Date 

UT 

Slit PA 

Exp. time 


dd/mm/yy 

hh:mm:ss 

[° E of N] 

M 

HD 97048 

21/03/05 

04:41:47 

0 

300 


21/03/05 

04:47:43 

0 

300 


22/03/05 

02:26:21 

0 

300 


22/03/05 

02:32:18 

0 

300 


27/03/05 

00:40:24 

0 

450 


27/03/05 

00:48:49 

0 

450 


27/03/05 

01:00:48 

45 

450 


27/03/05 

01:09:12 

45 

450 


27/03/05 

01:18:56 

90 

450 


27/03/05 

01:27:22 

90 

450 


27/03/05 

01:36:49 

135 

450 


27/03/05 

01:45:14 

135 

450 

HD 100546 

27/03/05 

01:58:59 

150 

105 


27/03/05 

02:01:40 

150 

105 


27/03/05 

02:05:25 

15 

105 


27/03/05 

02:08:06 

15 

105 


27/03/05 

02:12:03 

105 

105 


27/03/05 

02:14:44 

105 

105 


27/03/05 

02:18:36 

60 

105 


27/03/05 

02:21:16 

60 

105 


average the continuum over 5 spectral pixels. Hence we ob¬ 
tain an average spatial profile of the continuum blue- and red- 
ward of the [Oi] line. To determine the interjacent continuum, 
a linear interpolation in the spectral direction is done at each 
spatial pixel row. We end up with a 2D continuum. The nor¬ 
malized spectrum of the source at each sampled spatial posi¬ 
tion is then retrieved by dividing the flux in each spatial row of 
the 2D echelle frame by the corresponding row of the 2D con¬ 
tinuum. This procedure is the straightforward 2D translation 
of the “standard” normalizing procedure. Furthermore, a 2D 
continuum-subtracted spectrum is computed by subtracting the 
2D continuum from the 2D airglow-corrected science frame. 
In Fig. [I] we schematically summarize the airglow-removal and 
normalizing method described above. 

3. Analysis 

AVD05 have shown that the double-peaked profiles observed 
in the [O i] spectra of many Herbig Ae/Be stars emanate from 
excited oxygen atoms in the surface layers of an inclined, ro¬ 
tating flared disk. HD 97048 and HD 100546 are two rela¬ 
tively nearby (distances 103 and 180 pc for HD_100546 and 
HD 97048 respectively, Ivan den Ancker et~ . mi Herbig 
stars that were shown to display these double-peaked [Oi] 
emission lines in this study. If the emission region is spatially 
extended, it is possible to extract spatial information on these 
objects in the [Oi] line. When altering the position angle of 
the slit on the sky, different spatially dependent spectra are ex¬ 
pected to be observed. If the slit is positioned over the target 
along the major axis of the disk, the normalized spectrum on 
one spatial side of the order will primarily show the blue peak 
of the profile, while the spectrum on the other spatial side will 
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spectral pixels 



spectral pixels 



Fig. 1. The sky-subtracted 2D echelle frame ( top left panel) 
contains the spatially (vertical) and spectrally (horizontal di¬ 
rection) dispersed spectrum of the target. The logarithmic color 
scale is the same in the four panels and represents the mea¬ 
sured flux in each pixel. The bright horizontal band (order) is 
the observed flux of the target, while the narrow vertical band, 
which fills the entire slit, is the airglow of Earth’s atmosphere. 
The [O i] feature is located approximately in the middle of the 
2D frame. The two “dips” in intensity left and right of the 
feature are two telluric absorption lines. For all our spectra, 
these absorption lines are located outside the spectral range of 
[O i] feature. Therefore there is no need to correct for this at¬ 
mospheric effect. Based on the spatial rows “above” and “be¬ 
low” the echelle order (horizontal lines), the average spectral 
shape of the airglow feature can be determined. This average 
airglow feature ( bottom left panel) is subtracted from the 2D 
science frame. The spatial profile of the continuum blue- and 
redward of the [O i] feature (vertical lines) is used to compute 
the 2D continuum ( top right panel ) by linearly interpolating in 
spectral direction for each spatial row. Finally, the continuum- 
subtracted, airglow-corrected 2D science frame, which only 
contains the [Oi] feature, is extracted ( bottom right panel). 

The normalized 2D spectrum, which is not shown in this figure, 
is retrieved by dividing the airglow-corrected science frame by 
the 2D continuum. 


predominantly display the red peak. In the case of the slit ori¬ 
ented perpendicular to the disk’s major axis and centered on the 
star, no change in spectral shape at different spatial positions is 
expected: the effect is wiped out due to spatial integration in the 
direction perpendicular to the slit. Intermediate PAs of the slit 
on the sky lead to intermediate changes in spectral profiles. No 
effect of slit PA is expected outside the emission line, because 
the continuum consists of photospheric emission of the central 
star. 

Apart from the normalized spectra at different spatial posi¬ 
tions along the slit, we have deduced a second quantity which 
may help to interpret the observations. As a function of wave¬ 
length (or equivalently radial velocity) we have determined the 
spatial peak position of the echelle order in the continuum- 
subtracted 2D frame. The quantity is deduced by fitting a 


Gaussian function to the spatial profile in each spectral column 
(this technique is called spectro-astrometry. e.g. iTakami et al.l 
120041 and references therein). We have compared these values 
to the peak position of the underlying continuum. If the slit is 
positioned along the major axis, one expects to see the spa¬ 
tial peak position shift from one side of the continuum peak 
position to the other when running through the (spectral) line 
profile. Together with the spatial peak position, the Gaussian fit 
delivers a spatial full width at half maximum (FWHM). In the 
—spatially unresolved— continuum this quantity is a measure 
for the atmospheric seeing. In the [Oi] 6300A line the FWHM 
is potentially larger, if the extent of the emission region is sig¬ 
nificant compared to the atmospheric broadening. As can be 
derived from Table [3] this is not the case in our observations. 
The FWHM ratios of the feature and the continuum are within 
the error bars all equal to unity. We will come back to this point 
in Sect.E] 

Table 2. The ratios of the measured FWHM in the feature and 
in the continuum, 1 -WH M/ ra ,/FWH M r „,„. All ratios are within 
the error bars equal to unity. Based on our model assumptions 
(Sect. EJ, the expected maximum value for the ratio is 1.08 
with the target at a distance of 103 pc and 1.03 at a distance 
of 180 pc. The slit PAs for the models indicated in this Table 
are relative to the modeled disk major axis (e.g., PA=0° means 
slit and major axis aligned). 


Object 

Slit PA 
[° E of N] 

FWHM f „, 

[arcsec] 

FWHM /ra ,/FWHM co „, 

HD 97048 

0 

1.00 

1.00 ±0.04 


45 

0.93 

1.00 ±0.04 


90 

0.98 

1.01 ± 0.04 


135 

1.13 

1.00 ±0.04 

HD 100546 

15 

1.32 

1.01 ± 0.04 


60 

0.95 

1.02 ± 0.03 


105 

0.98 

1.01 ± 0.03 


150 

1.05 

1.02 ± 0.03 

AVD05 180pc 

0 

1.00 

1.03 


45 

1.00 

1.02 


90 

1.00 

1.03 


135 

1.00 

1.02 

AVD05 103pc 

0 

1.00 

1.08 


45 

1.00 

1.06 


90 

1.00 

1.08 


135 

1.00 

1.06 


In Figs. |3 and 0 we have shown the (spatially integrated) 
spectra of HD 100546 and HD 97048 at each slit PA (upper 
panels). The small differences in spectral shape are likely due 
to the small slit width (0.3") and the extended targets, which 
implies that only a part of the disk falls inside the slit. The 
second row of figures contains the normalized spectra at dif¬ 
ferent spatial rows around the central spatial position of the 
order. The latter is the pixel row which contains the peak of the 
spatial profile of the continuum. Note that, due to limited see¬ 
ing, the spatial information is blurred out over different spatial 
rows. Nonetheless, when moving away from the central spatial 
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row, the wings may become more pronounced compared to the 
rest of the feature, if a small offset of the blue or red emission 
region on the sky relative to the central star is present. The bot¬ 
tom panels display the spatial peak position of the continuum- 
subtracted [O i] feature relative to the peak position of the con¬ 
tinuum. We have included this plot because it is insensitive to 
asymmetry in the spectral line profile: for each velocity, the de¬ 
termined central peak position is independent of the total flux at 
that spectral position. These plots give a more precise view on 
the spatial peak position of the blue and red wing of the feature 
than the normalized spectra, in which the relative peak flux of 
the wings influences the spectral shape at all spatial rows due 
to atmospheric blurring. 

The normalized spectra and spatial peak positions in Figs. [3 
and|3]show the signature of an emission region which is linked 
to the rotating circumstellar disk. For HD 100546, the closest of 
both targets, this fingerprint is clearest. With the slit positioned 
150° east of north, i.e. along the disk major axis, the normal¬ 
ized spectra at the NW side of the star show the red wing more 
pronounced than in the total spectrum, while the blue wing is 
clearly dominant on the SE side. The same behaviour is seen 
in the shift in spatial peak position in the feature relative to the 
continuum: the photocenter of the red wing lies NW, and that 
of the blue wing SE of the star. The same characteristics are 
seen in the spectrum at PA 105°, although less strong. Due to 
the differential angle of 45° between the slit and major axis PA, 
the spatial extent of the [O i] emission region translates only 
partially into a spectral effect. At slit PA 60°, no difference in 
spatial peak position between the blue and red wing is noted. 
The normalized spectra show no trend from SW to NE either. 
However, the entire [O i] feature seems to be shifted to the NE 
side of the star. We will come back to this in Sect. 0 At PA 
15° east of north, the spatial profile is reversed: the slit is turned 
135° on the sky and the north now appears on the top side of 
the order. The peak position of the red wing is clearly on the 
N side, confirming the conclusion derived from the spectra at 
slit PA 105 and 150°. The blue wing in the PA 15° spectrum on 
the other hand does not shift significantly to the S of the central 
star. In Sect. El we discuss this anomaly. 

The distance to HD 9 7048 is estimated to be 180 pc 
Ivan den Ancker et aT.l 1998l l. The extent of the [O i] emission 
region is therefore less clearly present in the data. Nonetheless 
the spatial peak position in the spectra with PA=0° suggests 
that the red peak is located northward of the central star, while 
the blue peak lies to the south. This behaviour is confirmed by 
the spatial peak positions in the PA= 135° spectrum: the blue 
wing appears on the SE side, the red wing on the NW side. The 
two spectra in between (at PA 45° and 90° east of north) do not 
display a clear spatial effect. 


4. Discussion: confronting the AVD05 model 

4.1. The model as starting point 

In AVD05 we have overplotted the spectral profiles of the [O i] 
lines of HD 100546 and HD 97048 with our model results 
(Fig. 26 in that paper). In this Section we go back to these 
models and confront not only their spectral profile, but also 


their spatial extent with the UVES observations. The param¬ 
eters of the model which we consider here include a stellar 
mass of 2.4 M 0 , inclination i = 45° and a surface-density profile 
oc R 2 5 . 80% of the [O i] emission in our model emanates from 
the region between 0.8 and 20 AU. For details on the modeling, 
see AVD05. 

To simulate the effects of atmospheric blurring and instru¬ 
mental spectral broadening of UVES, we have convolved the 
theoretical (spatially extended) emission line with Gaussian 
profiles. In spatial direction, we have taken a FWHM equal to 
1"—a typical value derived from the UVES observations, see 
Table |2]— while in spectral direction we take the instrumental 
width to be equal to the average width of six (unresolved) tel¬ 
luric absorption lines near the [Oi] feature (3.7 km s 1 ). The 
results are presented in Fig. El The four top-panel spectra are 
identical, since we did not simulate the small slit width, but 
just integrated over the entire disk. The middle panels display 
the normalized spectra at each spectral row. The shape of the 
spectrum changes depending on the relative position angle be¬ 
tween slit and disk major axis (APA) and the spatial distance to 
the center of the target. The bottom panels show the theoretical 
peak position throughout the feature. 

Focusing on the bottom row of panels, the following re¬ 
marks can be made. Qualitatively, the effects seen in the mod¬ 
eled line profile agree with the observations, indicating that 
the rotation of the disk indeed causes the double-peaked [O i] 
6300A lines in HD 97048 and HD 100546. Nevertheless, there 
is a anomaly in the magnitude of the spatial effect: the shift 
of the peak positions in the modeled profile is approximately 
three times larger than the observed shift. The presented model 
is however not a fit to our data, but an ad hoc solution. As we 
have shown in AVD05, this solution is far from unique in ex¬ 
plaining the observed (spectral) profiles. Moreover, the effects 
of atmospheric and instrumental broadening have been mim¬ 
icked by convolving the profile with Gaussian curves. This is a 
first-order approximation. The real point spread function may 
blur out the shift in spatial peak position even more or at least 
differently. 

To show that the shift in spatial peak position in both targets 
is in qualitative agreement with the models, we have calculated 
the relative difference in peak position (DPP) between the red 
and blue wing. The spectral position of the wings is indicated 
with the vertical lines in Figs. d Eland 0 The relative peak po¬ 
sition is computed for each slit PA. For the models the values 
are normalized to unity by dividing them by the maximum DPP 
(at APA = 0°). For HD 97048 and HD 100546, we have fitted 
the DPPs with a cosine function DPP = A cosfPA - B ) with 
A and B fit parameters. The best-fit (y 2 , d ) amplitude A is used 
to normalize the DPPs to unity. The best-fit PA offset B indi¬ 
cates the position angle of the major axis of the disk. In Fig. El 
the modeled and observed DPPs are displayed. From the fit, 
we can derive that the PAs of the disk major axis in HD 97048 
and HD 100546 are respectively 160+19° and 150+11° east of 
north 2 . Note that the —independently determined— PA value 
for HD 100546’s disk is in perfect agreement with the known 

2 As a convention we have defined the PA of the disk major axis in 
the direction of the blueshifted part of the disk. 
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Fig. 2. Top panels: The spatially integrated spectra of HD 100546 at different slit position angles (PA=15°, 60°, 105° and 
150° east of north). The normalized flux is plotted versus the (relative) radial velocity. 

Middle panels: The normalized spectra at each spatial pixel row around the central spatial position. The spectra are shifted for 
clarity. The formal spatial distance between two adjacent spectra is 0.182", but atmospheric seeing blurs the spectra in spatial 
direction. 

Lower panels: The spatial peak position of the continuum-subtracted feature at each spectral row are plotted. The offset is relative 
to the spatial peak position of the underlying continuum. The scale is in arcsec. The curves are independent of the peak flux in 
each spectral row and hence give a clear view on the spatial position of each (spectral) part of the line profile. In the top and 
bottom right of the figures, the positions on the sky —in the direction of the slit PA— are indicated. 

The vertical dotted lines indicate the spectral position of the blue and red wings. 


PA of this system (^150° E of N, e.g. iLiu et al.l 120031 . We 
conclude that the S-SE parts of the disks around HD 97048 
and HD 100546 rotate toward us, while the N-NW side moves 
away from us. 

Note that the spatial peak position of the feature in 
HD 100546 is not symmetric. In the [Oi] line, the NE side of 
the order seems to be brighter than the SW side, pulling the 
peak position —which is a measure for the photocenter— in 
this direction. This asymmetry may be due to the fact that the 
disk around HD 100546 is a flared disk: the surface of the side 
of the disk closest to the observer is seen under a more acute 
angle than the surface of the opposite side. Therefore, the pro¬ 
jected surface of the nearest side of the disk is smaller than the 
projected surface of the opposite side. The modeling of this ef¬ 
fect is beyond the scope of the present paper. However, to give 


an idea about its magnitude, we can make a simple estimate. 
Assuming a typical, homogeneous flared disk with an open¬ 
ing angle of 20° and an inclination of 50°, the viewing angle 
between the line of sight and the normal to the surface of the 
nearest side of the disk is about 70° (i.e. almost “edge-on”). 
The viewing angle on the opposite side is approximately 30°. 
The projected surface of the nearest side of the disk is under 
these assumptions only ~30% of the projected surface of the 
other side of the disk. As opposed to the spatial shift of the 
blue- and redshifted side of the disk, which is maximal when 
slit and disk major axis are aligned, this effect influences the 
spatial peak positions in a spectrum most when the slit and ma¬ 
jor axis are perpendicular. When slit and major axis are aligned, 
the spatial integral —perpendicular to the slit— wipes out the 
signature of the differences in projected surface. This is consis- 
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PA = 0° PA = 45° PA = 90° PA = 135° 



tent with the observations of HD 100546: at PA 15° and 60°, 
the effect may play a role. This would indicate that the NE side 
of the disk is further away from the observer than the SW side. 
If so, our knowledge on the rotation of the [O i] emission re¬ 
gion can be further refined: not only can we deduce that the 
blue wing is located on the SE side of the disk, but also that the 
rotation occurs counterclockwise around the central star. 

A second effect which is not taken into account in our 
model is possible inhomogeneity in the [Oi] emission re¬ 
gion. Obviously, deviations from homogeneity do occur: the 
spectral [Oi] feature is not symmetric and temporal changes 
in the line profiles are observed. For both targets, we have 
other medium-to-high resolution 6300A spectra at our disposal 
(ESO 3.6m CES and ESO 1.5m FEROS, and CAT CES for 
HD 100546 only). These spectra are described in detail in 
AVD05. In Fig. |6] the spectra are displayed. The plot allow us 
to take a first glimpse at the temporal behaviour of the [O i] line 
in Herbig Ae/Be disks. The CAT spectrum (A/AA = 65,000) is 
taken in 1994, the FEROS spectra (4/Ad = 45,000) in 1999 
and the ESO 3.6m spectra {A/AA = 120,000) in 2002 (see 
AVD05 for the exact dates). Assuming that emitting gaseous 
“clumps” survive in the disk surface over multiple orbital pe¬ 
riods, the temporal changes in the spectral shape of the [O i] 
line must display periodic behaviour, linked to the disk’s ro¬ 


tation. The respectively 3 and 4 spectra of HD 97048 and 
HD 100546 are hardly sufficient to describe potential period¬ 
icity in the [Oi] 6300A line, but some conclusions can nev¬ 
ertheless be drawn. First of all, the equivalent width ( EW ) 
of the [Oi] 6300A line changes with time. For HD 97048, 
this change is small (|£W|=0.141, 0.135 and 0.160A in the 
FEROS, ESO 3.6m and UVES spectrum) while it is larger for 
HD 100546 (I.EWI =0.183,0.180,0.152 and 0.16lA in the CAT, 
FEROS, ESO 3.6m and UVES spectrum). Also the variations 
in spectral shape are larger in HD 100546. Due to the short 
time base of the observations —6 and 11 yr for HD 97048 and 
HD 100546 respectively— only the influence of fast orbiting 
clumps can potentially be detected: for a 2.4 M 0 star, 6 and 
11 yr correspond to the orbital period at 4.5 and 7 AU respec¬ 
tively. Inhomogeneities closer in orbit the star at least once in 
this time span. Due to the fact that spectroscopic observations 
are only sensitive to the projected velocity, non-homogeneous 
emission from the high-velocity part of the disk influences the 
spectral shape of the feature over the entire velocity range, 
and not only at high projected velocities. For HD 100546, the 
double-peaked profile seems to alternate between a blue-wing 
dominated (CAT, UVES) and a red-wing dominated feature 
(FEROS, ESO 3.6m). We interpret this as additional evidence 
for an [0 1 ] emission region which corotates with the circum- 
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APA = 0° APA = 45° APA = 90° APA = 135° 



projected velocity [km/s] 

Fig. 4. Similar plot as Figs.|2]and|3for the model. If the PA of the slit is aligned with the PA of the major axis of the inclined disk 
(APA = 0°), the spatial effect is most pronounced. When slit and disk major axis are perpendicular (APA = 90°), no spatial effect 
is seen. The parameters of the model presented in this figure consist of the stellar parameters of an AOV star, a surface-density 
power law 2 oc R 2 \ an inclination i = 45° and a distance to the source of 180 pc (see AVD05 for details on the modeling). 
The line profile and continuum have been broadened with a Gaussian curve in spatial (FWHM = 1") and spectral direction 
(FWHM = 3.7 km s 1 ) to simulate the dispersion due to the atmosphere and the UVES spectrograph. 


stellar disk. Furthermore, the presence of clumpy emission pro¬ 
vides an alternative explanation for the observed spatial peak 
positions at different slit PAs. We come back to this point in 
Sect. |5] when we discuss the possible origin of the clumps. 

4.2. The observations as starting point 

The accuracy of our model can also be tested starting from 
the observations. The spectral profile of the [Or] 6300A line 
contains information on the velocity field of the [Oi] emis¬ 
sion region. Assuming a Keplerian velocity field, a disk in¬ 
clination and a distance to the target, the spatial extent of 
the [Oi] emission region can be calculated from the emis¬ 
sion line. First we have averaged the blue- and red wing of 
the spectral profile in order to construct a symmetric 6300A 
line. The highest-velocity (v max ) bin which contains [O i] flux 
(i.e. where the wing starts to stand out above the continuum) 
can be linked to the innermost part of the [O i] emission re¬ 
gion. The observed flux in this bin is attributed to a ring at 


radius R = GM*(sin i/v max ) 2 . Only this inner ring contributes 
to the flux in the v' max bin, because the Keplerian velocity is 
smaller at larger radii and because no [O i] emission is detected 
at higher velocities. We subtract the theoretical emission profile 
of this ring, which covers all projected velocities from -v' max to 
+v max , from the total spectral line profile. The spectral bin at 
''max of the residual is now empty and we continue with the 
next high-velocity bin. The method is repeated until the spec¬ 
tral profile is completely dismantled down to v = 0 km s 1 . 
In this way the observed intensity-versus-radius profile can be 
reconstructed. For both FID 97048 and FID 100546 we have 
taken M* = 2.4 M e and i = 45°. In Fig. 0 we display the ob¬ 
served intensity-versus-radius curves, I(R), deduced from the 
spectral data of FID 97048 and FID 100546, and our modeled 
intensity distribution. We have multiplied the curves with the 
surface S (R) of the ring at radius R. In this way, the relative 
contribution at each radius can readily be compared. Since the 
modeled flared disk is truncated at the dust destruction radius 
(0.8 AU in the case of an AOV star), in our models the [Oi] 
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Fig. 5. The relative difference in peak position (DPP) between 
the blue and red wing for the model (solid line) and observa¬ 
tions (squares and diamonds) versus the relative position angle 
(APA) of slit and disk major axis. The observed DPPs are nor¬ 
malized by dividing with the amplitude of the best-fit cosine. 
A relative shift of 0 indicates that the spatial peaks of the blue 
and red wing of the feature are located at the same position. The 
best-fit (^" ed = 1.9) PA of the disk major axis in HD 97048 is 
160+19° east of north. For HD 100546, the best-fit (*f ed = 0.6) 
PA angle is 150+11° east of north, which coincides perfectly 
with the known PA of the disk. 


HD 97048 HD 100546 



Fig. 6. Temporal changes in the spectra of HD 97048 and 
HD 100546. Top panels: the available normalized spectra, 
shifted according to their relative age (0.1 corresponds to 
5 yr). From top to bottom: UVES, ESO 3.6m, FEROS and, for 
HD 100546 only, CAT. Bottom panels: the flux difference of 
the normalized UVES spectrum and the other normalized spec¬ 
tra, again shifted according to relative age. These differential 
spectra are rebinned to the resolution of the lowest-resolution 
spectrum of the two. From top to bottom: UVES - ESO 3.6m, 
UVES - FEROS and, for HD 100546 only, UVES - CAT. 


emission region also starts at this distance to the star. The outer 


limit is chosen to be 100 AU. Beyond this limit, the modeled in¬ 
tensity drops quickly and is insignificant compared to the [O i] 
intensity closer to the star. The observed I(R) x S (R) curve of 
HD 97048 has one prominent peak. Under the current assump¬ 
tions of stellar mass and inclination, 80% of the total [O i] emis¬ 
sion emanates from the region between 0.9 and 12 AU. The 
curve of HD 100546 is qualitatively equal to that of HD 97048 
out to 5 AU, but shows a strong second “bump” further out, 
around 15 AU. 80% of the total [Oi] emission emanates from 
the region between 1.2 and 23 AU, with 25% emanating from 
the bump. The observations suggest that the [Oi] emission re¬ 
gions of HD 97048 and HD 100546 are more concentrated 
around the central star than in our model. The drop-off at large 
radii is significantly stronger than predicted in both sources. 
The inner edge of the model (0.8 AU) on the other hand is in 
agreement with the observed intensity-versus-radius curve: no 
significant [Oi] flux is detected at high velocities. 

We have confronted the observational intensity-versus- 
radius curves derived from the spectral line profiles to the spa¬ 
tial extent deduced from the spatial information contained in 
the UVES spectra. Assuming a stellar mass of 2.4 M 0 , an in¬ 
clination of 45° and a distance of 180 pc for HD 97048, the 
spatial peak positions throughout the [Oi] line derived from 
the spatial and spectral information agree. The maximum value 
of the ratio 1 ; W H M f eat /l -WHM r „„, derived from the observed 
intensity-versus-radius curve is 1.01. This value is in agreement 
with the values in Tabled the spatial extent is not resolved in 
the spatial FWHM of the feature. Adopting the same stellar 
mass and inclination for HD 100546 and a distance of 103 pc 
leads to a significant difference in spatial peak position as de¬ 
rived from the spectrum-based intensity curve and the spatial 
UVES information. The spatial peak position estimated from 
the intensity curve is 1.5 times larger than the measured UVES 
peak position. The ratio of the FWHM in the feature and in the 
continuum is 1.05. This is on the verge of the spatial resolution 
of our UVES observations (see TableQ. This factor may arise 
from an erroneous distance to the source. The latter is needed 
to convert the radial scale to the angular scale on the sky. The 
possibility that the distance to HD 100546 is underestimated by 
a factor of 1.5 seems how ever unlikely, considering th e small 
error bars (±7 pc) given bv lvan den Ancker et al.l l 1 9981) . In or¬ 
der to reconcile both measurements, either the stellar mass of 
HD 100546 has to be increased to 3.6 M 0 or the inclination 
must be increased to 60°, or a combination of both proportional 
to M* sin 2 i. When incorporating the 1.5 factor, the maximum 
ratio of the FWHM in feature and continuum decreases to 1.03, 
undetectable in the current UVES observations. The first peak 
in I(R) x S (R) would then shift from 1.5-5 AU to 1-3 AU and 
the second [Oi] emission “bump” at 10-25 AU to 7-17 AU. 

5. A massive planet orbiting HD 100546? 

5.1. The presence of a gap reconfirmed 

iBouwman et alJ J2003I) have suggested that the mid-IR contin¬ 
uum emission in the SED of HD 100546 points to the pres¬ 
ence of a dust-poor cavity around 10 AU. This cavity is be¬ 
lieved to be planet-induced. Due to lack of dust at 10 AU, 
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the far-side rim of the gap is directly irradiated by the cen¬ 
tral star. The mid-plane temperature is therefore high and the 
vertical scale height of the disk at the far side of the gap is 
more than a factor of two that of a standard flared disk at the 
same radius. The suggested gap and wall explain the relative 
lack of mid-IR continuum flux shortward of 10 pm and the 
strong excess flux at longer wavelengths. The presence of this 
gap has been confirmed by mid-IR nulling interferometric mea¬ 
surements dLiu et al . 200* and Hubble Space Telescope (HST) 
STIS observations dCiradv et al.l2005l) . The latter authors have 
observed HD 100546 in the far-UV wavelength range, obtain¬ 
ing long-slit spectra. They suggest that the gap has an outer 
radius of 13 AU and is not centered on the star but rather on a 
region 5 AU southeast of the star. They exclude that an earlier- 
than-M5-type stellar companion is present within 2.5-10 AU of 
HD 100546. 

The intensity-versus-radius profile derived from the [Oi] 
feature in HD 100546 (Fig.0 may also be interpreted in this 
framework. The I(R) x S (R) curve in HD 100546 a local max¬ 
imum at short radii and declines furtherout up to ~5 AU. The 
curve of HD 97048 displays the same behaviour in this region, 
which suggests that the inner parts of the [O i] emission region 
are similar in both sources. Further out however, the disk struc¬ 
ture of HD 100546 appears to be significantly different than 
that o f other Herbig A e /Be st ars s uch as HD 97048. E xactly 
where iBouwman et al.l d2003l) and IClradv et all d2005l) locate 
the outer edge of the gap, there is a significant excess in [O i] 
emission in HD 100546. Note moreover that the I(R) x S(R) 
curve of HD 100546 is continuous and has a local minimum 
at 6-10 AU. The [Oi] emission region is suggested to be lo¬ 
cated in the disk surface layer (AVD05). If an additional “wall” 
(apart from the inner rim at 0.8 AU) is present in the disk, the 
irradiated disk surface is significantly larger than in the stan¬ 
dard flared-disk geometry. More [O i] emission than in a stan¬ 
dard flared disk would be observed. Although the dust disk may 
be mostly cleared at 10 A U, the gas disk is very likely still 


present around this radius feouwman et al.lEoQ3t iGradv et al.l 


2005). Therefore, the continuity of the [O i] intensity-versus- 


radius curve is not in disagreement with the possible presence 
of a dust-free gap in the disk. The local minimum in the in¬ 
tensity curve is located at the radius where the “standard” [O i] 
emission region ends and the “additional” wall emission be¬ 
gins. 

The gap in the disk around HD 100546 is likely planet- 
induced. The mass of the planet can in principle be deduced 
from the size of the gap, the planet’s orbital radius and the stel¬ 
lar mass: according to lArtvmowiczi ( 1987 1 the diameter of the 
gap cleared by a planet is equal to 4 V3 Rr^cUc, where RR 0C he is 
the radius of the planet’s Roche lobe. An additional constraint 
on the planetary mass can be obtained fro m il .in & PanaloizoiJ 
J 1 9931) . The authors state that a g ap can only be created if the 
radius of the planet’s Hill sphere dHi I III 18781) is larger than the 
vertic al scale height of t he disk . In their estimate of the planet’s 
mass, IBouwman et al.l d2003h assum e a relative scale height 
H/R = 0.045 at 10 AU, b ased on thelDullemon det^p J20()ll) 
model for AB Aur. In our lChiang & Goldreichl ill997ft flared- 
disk mod el for HD 100546, H/R is equal to 0.071 at 10 AU 
(AV D05 ). IBou wma 11 et al. 12003) double the H/R value argu¬ 



1 10 100 

R [AU] 



Fig. 7. The intensity-versus-radius curve I(R) derived from the 
spectral shape of the [Oi] 6300A line in HD 97048 (top) and 
HD 100546 (bottom), and the model curve. The intensities have 
been multiplied by the surface S (R) of the ring at radius R. The 
fluxes of the inner and outer parts of the disk can be directly 
compared in this way: the total strength of the feature is the 
sum of the plotted values. The error bars have been deduced 
from the signal-to-noise ratio of the spectra. The model curve is 
only non-zero in the region between the dust destruction radius 
at 0.8 AU and the outer radius of 100 AU considered in the 
model. The observations suggest that the [0 1 ] emission regions 
in HD 97048 and HD 100546 are more compact than the model 
suggests. The peak of the I(R) x S (R) curve in HD 97048 is 
located between 1 and 5 AU. For HD 100546 a second peak 
can be discerned at 10-25 AU. See text for explanation. 


ing that the vertical scale height of the wall at 10 AU must be 
twice as large as predicted by the flared-disk models in order to 
explain the mid-IR continuum fluxes. However, the suggested 
wall develops due to the direct radiation of the central star. The 
stellar photons can only reach the edge on the far side of the 
gap if a gap is already present. In our opinion it is therefore 
not necessary to assume a larger vertical scale height of the 
disk than in the original flared-disk model. The lower limit to 
the companion’s mass is then M p i anet > (H/R) 3 3 M* = 0.7 or 
2.7 Mj up i tel depending on the chosen flared-disk model. 
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Fig. 8. The size of the planet-induced gap in the disk of 
HD 100546, 4 V3 R Ro rhp, as a func tion of the planet’s mass 
jArtvmowiczl 1 1 987t Eggletonlll983ll . The curve represents a 
companion at a 6.5 AU distance from the central star. The 
horizontal dashed line is an estimates for the diameter of the 
gap (4.3 AU). The curves intersect at a companion’s mass of 
20 Mjupiter. The left vertical dashed line and arrow indicate the 
lower limit to the planet’s mass, 2.7 Mj up i ter , derived from the 
requirement that the planet’s Hill sphere radius is larger than 
the local vertical scale height in the disk. The vertical line on 
the right side indicates the total disk mass in small dust grains, 
30 Miiinifp,-, derived from the 1.3mm flux (e.g. lAcke et all20 04). 
iGradv et all J2005I) exclude th e presence of an earlier -than-M5 
type companion. According to Schmi dt-Kaleil (1 19821) . this cor¬ 
responds to an upper limit on the mass of ~ 30 Mj up i ter . We 
have indicated this upper limit with the second arrow. 

In Fig. [x] the expected diameter of the gap is shown as a 
function of the planet’s mass. The curve is based on the ap- 
proximation for the Roche lobe radius described by tEggletoiJ 
( 1983 ). The radius of the planet’s orbit is assumed to be the 
location of the local minimum in the I(R) x S (R) curve of 
HD 100546, about 10 AU. We have estimated the size of the 
gap by determining the FWHM around the local minimum. 
The gap diameter we derive is 6.5 AU. Comparing the the¬ 
oretical and observed gap diameter, a mass of 20 Mj llp j te r is 
deduced for the companion. An object with a mass exceed¬ 
ing s s 13 Minnitpr lies beyond the deuterium-burning thresh¬ 
old lOpnenheimer et al.l 12000 ) and thus the companion of 
HD 100546 is a brown dwarf. Note that the spatial extent of 
the [Oi] emission region in HD 100546 estimated from the 
HR) x S (R) curve —based on the spectral information in the 
UVES data— differs by a factor 1.5 from the extent derived 
from the spatial information in the UVES data. If we apply this 
factor, the distance between planet and star becomes 6.5 AU 
and the gap diameter 4.3 AU. The estimated mass remains 
20 M Jup i ter however, since also R Roc h e reduces by a factor of 
1.5 in this case. An upper limit to the companion’s mass c an be 
derived from the observational constraints set bv iGradv et al.1 
J2005I) . They exclude the possibil ity of an earlier-than- M5 type 
stellar companion. According to Sch midt-Kaledll 19821) . an M4- 


type star has a mass of 0.3 M Q or equivalently 30 Mj up i ter . We 
have indicated this upper limit in Fig. 0 As a reference, the 
total disk mass derived from the 1.3mm flux of HD 100546 is 
also 30 Mj up i ter dAcke et all2004l) . 

5.2. Implications of the presence of a companion 

The rough estimates derived in the previous paragraph shed a 
new light on other measurements of HD 100546. A massive 
companion disrupts the original density distribution of the cir- 
cumstellar disk, stirring up the surrounding matter. Spots of 
over- and underdensities are bound to form, their orbits syn¬ 
chronized with that of the companion. The [Oi] emission re¬ 
gion in the environment of the companion is therefore expected 
to deviate from axisymmetry. The temporal changes in the line 
profile of [Oi] 6300A line shown in Fig. [6] appear to display 
variability with a typical timescale of ~ 11 yr. The suggested 
giant planet at 6.5 AU orbits the central star in close to 11 yr. 
The “clumps” in the [0 1 ] emission region of HD 100546 men¬ 
tioned in Sect. n~n may thus very wel l be in resonance with 
the orbiting planet. iGradv et all ( 2005 ) derived from the HST 
STIS long-slit spectra that the cavity is centered 5 AU southeast 
of the star. The spectra were obtained on September 22, 2001 
and June 23, 2002, i.e. approximately 6 months before and 
3 months after our ESO 3.6m CES spectrum of HD 100546. 
The latter spectrum shows a [0 1 ] 6300A profile with a domi¬ 
nant red wing. Combining this information with the newly ob¬ 
tained knowledge on the disk’s rotation, most of the [0 1 ] emis¬ 
sion emanated from the northwest part of the disk at that time. 
Assuming that the cavity is centered on the giant planet’s posi¬ 
tion hence leads to the conclusion that the gaseous clump which 
causes the temporal changes in the [0 1 ] profile is located on the 
opposite side of the star as seen from the companion. 

From the shape of the [0 1 ] 6300A line in the 2005 UVES 
spectrum of HD 100546, the clump must mainly be located in 
the blueshifted part of the disk. Nevertheless, the blue wing is 
not so dominant as the red wing in the 2002 ESO 3.6m spec¬ 
trum. This may indicate that the clump has not reached its high¬ 
est projected velocity or equivalently has not yet arrived on the 
disk’s major axis at PA=150° east of north. This image is rein¬ 
forced by the shape of the spatial peak position through the fea¬ 
ture at different slit PAs (see bottom panels Fig. [2}. The spatial 
peak position throughout the profiles at PA=15° and 60° east 
of north appear to be oriented more to the northeast/east than 
expected, especially in the blue wing. The deviation of the spa¬ 
tial peak positions from the modeled peak positions at these 
PAs can be readily explained by the presence of a [0 1 ] emit¬ 
ting clump between position angles 60° and 105° east of north. 
The reasoning outlined above predicts that the planet is at 
present (i.e. 2005) to the southwest/west of the star, between 
PA=240° and 285° east of north. The spectral shape of the [0 1 ] 
6300A emission line thus provides a powerful tool to prognos¬ 
ticate the companion’s position on the sky. 

The fact that the [Oi] emitting clump in the disk of 
HD 100546 is at the time of the UVES observations (2005) to 
the northeast/east of the star, and that the spectral profile was 
red-wing dominated in the ESO 3.6m spectrum (2002) provides 
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further knowledge on the rotation of the [O i] emission region 
and thus of the disk. In three years time, i.e. approximately one 
quarter of the estimated orbital period, the clump has moved 
from the redshifted region of the disk (i.e. the north/northwest) 
to the north/northeast of the star. This indicates that the emis¬ 
sion region orbits the central star counterclockwise. Note that 
we have already reached a similar conclusion for the entire disk 
in Sect. rrn However, the argumentation in that Section was 
based on the assumption that the flaring of the disk’s geom¬ 
etry causes the deviation in the spatial peak positions of the 
UVES data. In the present Section we have suggested that the 
signature in spatial peak position can also be attributed to the 
presence of a clump in the [O i] emission. Even though both ar¬ 
guments may not be valid side by side, the resulting conclusion 
(i.e. the disk’s rotation occurs counterclockwise) is the same. 
We believe that the presence of an emitting clump provides the 
best explanation for the whole of phenomena which give the 
[Oi] 6300A line its (temporal, spectral and spatial) shape. The 
influence of the flaring of the disk is probably less important. 

The outer parts (R ~ 300 AU) of the disk around 
HD 100546 display spiral-arm structure (jciradv e t al J 1200 it) . 
In an effort to explain this structure. lOuillen et al.l ( -2005 ) have 
modeled the effects of a massive gravitationally bound object 
circling HD 100546. They conclude that the mass of such an 
object needs to be of the order of 20 Mj up j t er or larger to explain 
the observed impact on the disk structure. They argue that the 
presence of such a massive companion at R <300 AU is un¬ 
likely, because it would have been detected in the HST images. 
They conclude that a recent fly-by of a nearby star is prefer¬ 
entially the best explanation for the spiral arms, although they 
cannot find a good candidate in the field. It is striking however 
that our simple estimates for the companion’s mass and orbital 
radius agree with the boundary conditions set bv lOuillen et alJ 
( 2005 8): the mass of the companion is sufficiently high to cause 
the spiral pattern in the outer disk. Furthermore, the compan¬ 
ion’s orbit is so close to the central star that the residuals of 
the stellar point spread function (PSF) would easily cover up 
its presence in all existing observations (e.g. lAugereau et all 

EqqI)- 

The presence of the 20 Mj up i ter companion described in 
this Section leaves a signature in the radial velocities of the 
central star. The mass center of the star-planet system lies at 
a distance of 0.12 AU from the star. The amplitude of the 
stellar radial-velocity variations caused by the massive planet 
hence is 0.33 km s _1 x sin i = 0.25 km s _1 assuming i = 
50°. This amplitude is rather small, considering the observa¬ 
tional difficulties. The optical spectrum of HD 100546 dis¬ 
plays very few photospheric lines. Furthermore, these absorp¬ 
tion lines are broadene d due to th e stellar rotation (v sin i = 
55 km s~*. I Acke & Waelkensii20()4j) . Therefore such a compan¬ 
ion would easily go undetected in the existing spectroscopic 
data on HD 100546. 

With the current instrumentation available, it seems impos¬ 
sible to directly observe the companion of HD 100546. The 
orbital radius of 6.5 AU corresponds to an angular separa¬ 
tion of 65 mas at 103 pc. Space-based telescopes have diffrac¬ 
tion limits which do reach this kind of high-angular resolution. 
Adaptive Optics (AO) instruments on ground-based 10 m tele¬ 


scopes are able to obtain this spatial resolution. The main prob¬ 
lem however is the large brightness contrast between star and 
planet. Even in more favorable parts of the near-IR spectrum 
(e.g. in a methane band), this brightness ratio is out of reach for 
the present-day instruments. However, the companion should 
be readily detectable with high-contrast AO instruments on fu¬ 
ture ground-based Extremely Large Telescopes (ELTs). 


6. Conclusions and discussion 


In this paper we have provided further evidence that the [O i] 
6300A emission in the disks around HD 97048 and HD 100546 
emanates from the disk surface, as was first suggested in 
AVD05. Thanks to the combination of high spectral resolution 
and spatial information on the [O i] emission region contained 
in the UVES spectra, we are able to determine the the blue- 
and redshifted part of the disk in the two targets. Although the 
emission region was not resolved in the sense that the FWHM 
in the feature is larger than the continuum FWHM, the observa¬ 
tions clearly show a variation in spatial peak position in the fea¬ 
ture. The dependence of that variation on the slit position angle 
agrees qualitatively with our [O i] 6300A emitting flared-disk 
model. Based on the observed intensity-versus-radius curve of 
HD 97048 and HD 100546, we find that the [O i] emission re¬ 
gion in our model is more extended than the emission region in 
the targets. 

The major axis of the well-studied circumstellar dust disk 
around HD 100546 has a known position angle of 150° east 
of north. We have independently determined the same PA from 
our data. Note that two different diagnostics (the scattering dust 
particles and the oxygen gas responsible for the [O i] line) lead 
to the same angle. This indicates that the [Oi] emission region 
is linked to the geometrical distribution of the dust, in agree¬ 
ment with the conclusions of AVD05. For HD 97048 we find a 
major axis PA of 160° east of north, which is the first determi¬ 
nation of this angle in the literature. 

Furthermore, we have obtained knowledge on the rotational 
direction of both disks. The disk of HD 97048 is oriented ap¬ 
proximately north-south with the south part rotating towards us 
and the north part redshifted. The disk of HD 100546 is posi¬ 
tioned following a NW-SE line, with the SE part blueshifted 
and the NW part moving away. For the latter object we have 
shown that the rotation occurs counterclockwise around the 
central star. The data for HD 97048 do not allow us to make 
a similar statement. 

For HD 100546, we find evidence for the presence of a gap 
in the circu mstellar disk at about 10 AU, as was first suggested 
bvlBouwman et all d2003l) and later confirmed bv lGradv et alJ 
(2005). The radial intensity distribution of the [Oi] emission 
region indicates the presence of excess emission at R > 10 AU 
in comparison to HD 97048. We suggest that this additional 
emission emanates from the region on the far side of the gap 
in the dust disk, where a wall has formed. The gap is likely in¬ 
duced by a massive planet. Based on simple assumptions, we 
derive a diameter for the gap of - 4 AU and an orbital radius 
of 6.5 AU for the planet. The mass of this object is estimated 
to be 20 Mj up iter. The companion’s orbital period (~ 11 yr) is 
synchronized with the temporal changes observed in the [O i] 
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line profile. We suggest that the deviations from axisymmetry 
in HD 100546’s [Oi] emission region are also planet-induced. 
The [Oi] line profile variations in HD 97048, which displays 
no evidence for the presence of a massive planet, are far less 
pronounced. The derived mass and orbital elements of the com¬ 
panion of HD 100546 are in agreement with the conclusions of 
Ouillen et al. ( 2005 1. The object cannot be directly detected in 


HST images, and is massive enough to produce the large-scale 
spiral-arm structure observed in the outer disk of HD 100546. 
The estimates derived in the present paper can serve future ef¬ 
forts to detect the companion orbiting HD 100546 by direct 
imaging with ground-based Extremely Large Telescopes. 

The [O i] 6300A emission line in HD 100546 may prove to 
be a valuable tool to probe the companion’s orbital parameters. 
We will set up a long-term monitoring campaign on the Swiss 
Euler 1.2m telescope in La Silla (ESO) to study the variability 
of the line in more detail. The results of this study can provide 
further refinements in our knowledge of the giant planet’s or¬ 
bital period and hence distance to the star. Monitoring the tem¬ 
poral changes of the [O i] line may also lead to a first estimate 
of the eccentricity of the planet’s orbit. 

The observations of HD 100546 mentioned in the present 
paper all point to the presence of a massive object orbiting this 
target. Although the estimated mass of the companion suggests 
that it is a brown dwarf, it appears to have a formation his¬ 
tory very alike that of a planet. First, the object is expected 
to be close to HD 100546, which is surrounded by a mas¬ 
sive disk. Second, based on the mass estimate (20 Mj up i ter ) the 
companion would appear in the “brown dwarf desert”: a lack 
of binary companions in the mass range between 10 Mj up n el 
and late-M-type starsJ^ 1 gjMcCarth y_&_Zuckermanl2004l and 
references therein). lArmitage"& Bonnelli (120021) ' propose that 
brown-dwarf companions form contemporaneously with the 
primary. If the disk mass is sufficiently high compared to 
the brown dwarf mass, disk-companion interactions make the 
brown dwarf migrate inward. This occurs rather rapidly (~ 
1 x 10 5 yr for a massive disk). Given that the observed mass 
of the disk around HD 100546 is at present still ~ 0.3 M e 
fe.g. lAcke"etal.l2004 . the initial disk mass should have been 
high enough to ensure rapid inward migration of the compan¬ 
ion. Furthermore, there does not seem to be an underabundance 
of single stars throughout the brown dwarf mass range. Single 
stars form due to self-gravitational contraction. Given the dis¬ 
crepancy between the presence of single and long-lived com¬ 
panion brown dwarfs, is unlikely that brown dwarfs in disk- 
dominated environments are formed following the same mech¬ 
anism of self-gravitation. 

On the other hand, a 20 Mj up i ter object in a massive disk 
may have a similar formation history as a giant planet. Due 
to gravitational instabilities in the disk, the object can gain 
mass and grow (e.g. iBossl 12005 ). If this mechanism is in¬ 
deed at work in massive protoplanetary disks, it appears log¬ 
ical that it favors low-mass, planet-like objects over somewhat 
higher-mass brown dwarfs. Together with the fast removal of 
initially formed brown dwarf companions, this would explain 
the brown dwarf desert which is observed in multiple sys- 
tems. The HD 100546 syste m is at least 10 million years old 
Ivan den Ancker et alJl998l) . The presence of a massive object 


is therefore more logical under the assumption of a formation 
which occured in the disk (long) after the formation of the cen¬ 
tral star. We conclude that the object orbiting HD 100546 is a 
rare specimen of a brown dwarf companion with a formation 
history which is similar to that of a regular planet. 

The disk of HD 100546 is peculiar and extraordinary in 
many ways. The longevity of the disk is striking in compari¬ 
son with other pre-main-sequence stars. Most disks dissipate 
within a few 10 5 yr, while the disk around HD 100546 is still 
clearly present. Also the presence of relatively strong accretion 
is exceptional in such an old system (Balmer line e miss ion and 
accretion events d etected in the UV, e.g. IPeleuil et aiJ l2004: 
iGradv et~ai~ll200 5 ). The presence of a planet in the inner parts 
of the disk would cause a flow of gas which is “funneled” to¬ 
wards the central star. Furthermore, the solid-state features ob¬ 
served in the near-to-mid-IR suggest a highly processed dust 
component, quite similar to the spectra of comets in our solar 
system (see e.g. the Infrared Space Observatory data described 
by iMalfait et al.. 1998,) . The degree of crystallinity of the dust 
disk makes HD 100546 an outlier in the group of Herbig Ae/Be 
stars iAcke&Tvan den Anckeil5004 l. The grain size distribu¬ 
tion deduced from the (sub-)mm slope of the spectral energy 
distribution indicates the majority of the emitting dust grains 
have small sizes (<sl mm), while the dust particles in a lot of 
younger systems underwent substantial grain growth (>1 mm, 
lAcke et alJl2004j ). Now that it becomes more and more clear 
that a 20 Mj up it er companion is present, it is likely that the char¬ 
acteristics which make HD 100546 special all find their roots 
in the interactions between the disk and this companion. Such 
interactions clearly maintain the disk. Based on the presence of 
small and highly crystalline dust particles, we suggest however 
that the observed disk is not pristine, but rather more evolved. 
The interactions of the companion and the larger grains and 
planetesimals in the disk —which are expected to have formed 
in an earlier stage— could tear apart the latter. The small grains 
which emanate from such destructive interferences then replen¬ 
ish the small grain population in the disk. The occurance of 
events like these can provide a satisfying solution to the un¬ 
usual longevity of the disk in some young stellar systems. 

A final pecularity which is worth mentioning in this context 
concerns the unusual photospheric abundances of HD 100546. 
HD 100546 appears to be a so-called A Boo star: the photo¬ 
sphere displays a deficit in metals (Mg, Si, Fe and Cr), while 
nitrogen and oxygen ha ve approximately solar abundances 
dAcke & Waelken jl2004l). HD 100546 is the only Herbig star 
in the sample of iAcke & Waelkensi which clearly displays this 
behaviour. Selective accretion of metal-depleted gas is com¬ 
monly inferred to explain this selective depletion pattern. The 
metals condense into dust particles at much higher tempera¬ 
tures than the CNO elements. The radiation pressure on dust 
grains prevents the latter to fall in, while the gas-state elements 
are accreted. The gas and dust should however be decoupled 
in order to prevent the dust to be dragged along onto the star. 
The presence of a giant planet in the inner parts of the disk 
around HD 100546 could create the specific conditions needed 
to produce the selective accretion and consequent photospheric 
depletion pattern. If most metals are locked inside dust particles 
(or even planetesimals), the gas that is “funneled” by the planet 
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from the disk onto the central star is metal poor. Intriguingly, a 
deficit of metals in the stellar photosphere is opposite to what 
is typical of most giant planet systems found to date, in which 
the stars are primarily metal rich. 

Many observational characteristics of HD 100546 can be 
explained by the presence of a giant, nearby planet. If the argu¬ 
mentation in the present paper is confirmed by future observa¬ 
tions, the companion of HD 100546 is the youngest exo-planet 
discovered until today. The constraints deduced from a study 
of such a young planet+disk system will prove to be extremely 
valuable in both the modeling of disk-planet interactions, as 
well as in improving our comprehension of disk evolution and 
planet formation. 
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